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Abstract 

The nonlinear sigma model for which the field takes its values in the coset space 
0(l,2)/0(2) X Z2 is similar to quantum gravity in being perturbatively nonrenormal- 
izable and having a noncompact curved configuration space. It is therefore a good 
model for testing nonperturbative methods that may be useful in quantum gravity, 
especially methods based on lattice field theory. In this paper we develop the theo- 
retical framework necessary for recognizing and studying a consistent nonperturbative 
quantum field theory of the 0(l,2)/0(2) x Z2 model. We describe the action, the 
geometry of the configuration space, the conserved Noether currents, and the current 
algebra, and we construct a version of the Ward-Slavnov identity that makes it easy to 
switch from a given field to a nonlinearly related one. Renormalization of the model is 
defined via the effective action and via current algebra. The two definitions are shown 
to be equivalent. In a companion paper we develop a lattice formulation of the theory 
that is particularly well suited to the sigma model, and we report the results of Monte 
Carlo simulations of this lattice model. These simulations indicate that as the lattice 
cutoff is removed the theory becomes that of a pair of massless free fields. Because 
the geometry and symmetries of these fields differ from those of the original model we 
conclude that a continuum limit of the 0(l,2)/0(2) x Zi model which preserves these 
properties does not exist. 
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1 Introduction 



Nonlinear sigma models in 4-dimensional spacetime are not perturbatively renormalizable. 
They are therefore often viewed as being at best phenomenological, incapable of representing 
physical reality at a fundamental level. However, perturbative nonrenormalizability has never 
been proved to imply nonrenormalizability, i.e., the impossibility of giving precise meaning to 
a model, even in the extreme ultraviolet limit, in some nonpcrturbative sense. The purpose 
of this and the paper following is to present evidence that it is indeed impossible to give 
such meaning to at least one model: the 0(l,2)/0(2) x Z2 sigma model. 

Interest in this model stems from the following features that it shares with quantum 
gravity: 

1. Each theory has a single scale constant \x having the dimensions of mass. 

2. Each has a noncompact curved configuration space. 

3. Corresponding Feynman graphs for each (in perturbation theory) have identical degrees 
of divergence. 

4. Neither is perturbatively renormalizable. 

Despite the negative result, experience gained in studying the 0(l,2)/0(2) x Z2 model is 
useful in approaching the more difficult problem of quantum gravity, not only in regard to 
methodology, but also in setting ranges for certain parameters and in discovering the limits 
of present supercomputers. Although a positive result, i.e., demonstration of the existence 
of a continuum limit for the 0(l,2)/0(2) x Z2 model, would have been of obvious interest, 
no conclusion can yet be drawn about the existence or nonexistence of a continuum limit 
for quantum gravity. In no part of our analysis have we found reasons to have expected a 
negative result a priori. Each theory requires separate investigation, and it must be stressed 
that quantum gravity differs in several important respects from a sigma model: 

1. Its invariance group, the diffeomorphism group, is local, whereas the invariance group 
0(1, 2) of the sigma model is global. 

2. The ground state of the sigma model is degenerate, and spontaneous symmetry break- 
ing occurs. There is no evidence that the vacuum state of quantum gravity is degen- 
erate. 
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3. Quantization of the gravitational field smears the light cone. The light cone for the 
sigma model remains fixed and sharp. 

No method is presently known for studying nonlinear sigma models nonperturbatively in 
the extreme ultraviolet other than that of direct assault via the functional integral of Feyn- 
man. In field theory no nonperturbative definition of the functional integral is known other 
than that of taking the continuum limit of the corresponding integral on a lattice. In this pa- 
per we describe the continuum theory, pretending that it exists, in order to motivate our later 
choice of lattice action and of things to compute. The theoretical framework that will guide 
the computations is given in some detail because there are important differences between 
noncompact sigma models and the more familiar compact ones. For example, noncompact 
models undergo no transition to a phase of unbroken symmetry. In the following paper we 
introduce a lattice action that is particularly well suited to models with curved configuration 
spaces, and we describe the lattice Monte Carlo simulation of the 2)/0(2) x Z2 model, 
using this action. 

2 Action 

When we pass to the lattice simulation on the computer we shall need to work in Euclidean 
space, but for most of this paper we shall work in a Minkowski space with metric {j]^iv) — 
{jf^) — diag(— 1, 1, 1, 1). The action functional of a general sigma model has the form 



where commas followed by Greek indices denote differentiation with respect to the Minkowski 
coordinates a;'^(/i = 0, 1, 2, 3) and can be raised and lowered by means of the Minkowski 
metric. The 0* are scalar fields that take their values in charts of configuration space, and 
the Qij arc the components (in those charts) of the metric tensor of configuration space. 
The configuration space itself is always a coset space of a Lie group, and the metric tensor 
is the natural group invariant one. The configuration space is therefore a symmetric space 
with a covariantly constant curvature tensor Riju and a strictly constant curvature scalar 
R{= Rij^^). Only those coset spaces are chosen for which the natural metric is positive 
definite. 

The fields 0* will be taken dimensionless. In units with h = c = 1 the constant fi has 
then the dimensions of mass. It is sometimes convenient to absorb the constant /i into the 
fields by redefining 




(2.1) 
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and to follow the more standard convention of regarding scalar fields in four dimensions 
as having the dimensions of mass. This redefinition is useful, for example, in a chart in 
configuration space in which the metric takes the form 

g,, = - \R%i<\>''<\>' + ■ ■ ■ , (2.3) 
^ikji being the curvature tensor at = in this chart. The action ( |2.1] ) becomes 

^=-\\ [^n - + ■ ■ ■ ) i^^' d'x , (2.4) 

and the role of coupling constant in the theory is seen to be played by the inverse quan- 
tity /i'^. Of course, the Lagrangian of the actions ( p.lD and ( p.4|) is generally nonpolynomial, 
and there is in fact an infinity of coupling constants: the powers of We prefer to view 
fi as simply a scale parameter for the theory. In recognition of the model's superficial resem- 
blance to quantum gravity we shall call fi the bare Planck mass. 

Many years ago Palais and MostowQ showed that coset spaces of semisimple Lie groups 
can always be embedded in vector spaces in such a way that the group actions on each 
coset space can be represented as linear homogeneous transformations of the corresponding 
embedding space. This means that the field variables 0* can be replaced by a larger set 
of variables 0" together with a set of Lagrange multipliers A"^ that enforce the embedding 
constraints. The A"^ themselves remain invariant under the actions of the group. Since only 
a single chart is needed in the embedding space it is often convenient to present nonlinear 
sigma models in this form. 

The 0(l,2)/0(2) x Z2 model is defined, in this form, by the action 

S = -i/"' / [^ab r,>. 4>'r + A (jlab + 1)] d^x, a, h e {0, 1, 2} , (2.5) 

{r]ab) = diag(-l,l,l), (2.6) 
together with the constraint 

0° > 1. (2.7) 

The embedding space is seen to be a 3-dimensional Minkowski space on which the actions 
of 0(1,2) are Lorentz transformations. The Lagrange multiplier enforces the constraint 

r/afe0'^0' = -l (2.8) 



1 R. S. Palais, J. Math. Mech. 6, 673 (1957). G. D. Mostow, Annals of Math. 65, 432 (1957). 
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which, together with (pI7|), identifies the configuration space as the upper sheet of a spacehke 
hyperboloid in the embedding space. 

The group invariant metric on the hyperboloid is that induced by the metric rjab of the 
embedding space. If one parameterizes the configuration space by 

0° = cosh s , 

0^ = sinh s cos 9 , 

(j? = sinh s sin ^, (2.9) 

then one may write the induced fine element in the form 

- (c/0°)% [d(f)^y + (d(l?Y = ds^ + sinh^ sdO'^ , (2.10) 

and the action ( p.l|) becomes 

S = j (^^^^ ^ M + sinh^ s 9^^ 6^^) d^x . (2.11) 

The coordinates s, 9 are here allowed to cover the configuration space an infinite number of 
times so that the fields s, 9 may be taken differentiable wherever the 0" are differentiable. 



3 Geometry of the configuration space 

The configuration space of the 0(1, 2)/0(2) x Z2 model has constant negative curvature and 
is noncompact. In fact it is topologically IR^ and can itself be covered by a single chart, for 
example by the components of the 2- vector 

The constraints (2.7) and (|2.8|) can be expressed by 



l + 4>\ (3.2) 
and it is not difficult to verify that the metric in this chart takes the form 

0'0J' ^ ^ , , 

9ij = Sij - , , ,2 ^ hj e {1, 2} . (3.3) 

1 + 
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The curvature tensor is also easily derived: 

Rijkl = — {gik djl — 9il 9jk) , 
Rij Rikj 9ij 1 

R = = — 2 . 

Most nonlinear sigma models that have been studied have compact configuration spaces. 
We shall see later that noncompactness of its configuration space leads to profound differences 
in the behavior of the 0(1, 2)/0(2) x Z2 model from that of the compact models. This will 
refiect itself in the need to deal with the lattice simulation of the model by methods that 
will in part be unfamiliar. We note at this stage just one fact: The area of a large circle 
in the 0(l,2)/0(2) x Z2 configuration space increases not quadratically but exponentially 
with radius. This means that the functional integral must probe an increasingly huge region 
of the space of field histories as the continuum limit is approached. 

We shall need later an expression for the geodetic distance A, in the configuration space, 
between two points having coordinates (s, 6) and (s', 6') respectively. The simplest way to 
compute this is to note first that s is the geodetic distance from the point </> = to the point 
{s,6). Therefore, since the geodetic distance between two points remains invariant under 
the actions of 0(1,2), one may apply a Lorentz boost that brings one of the two points to 
= 0. This is most easily done in terms of the variables (f)"', a G {0, 1, 2}, and one finds 

A = cosh-'[-r]abr<l>") ■ (3.7) 
Substitution of expressions (2.9) into ( p.7| ) yields 

A = cosh^^ [cosh s cosh s' — sinh s sinh s' cos{9 — 9')] . (3.8) 



(3.4) 
(3.5) 
(3.6) 



4 Conserved currents 

The infinitesimal actions of the invariance group on the variables 0* may be expressed in the 
form 

S^' = Q\{(P)SC , (4.1) 

where the are infinitesimal group parameters and the Q^a are the components (in the 
configuration space chart) of a set of vector fields on the configuration space, which 
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satisfy the Lie bracket relations 

[Qa,Qp] = -QjC\p, (4.2) 

the c'a/s being the structure constants of the group. The statement that the Lagrangian 

L = -lfi'g,,(j)\,(f)^,>^ (4.3) 
is invariant under ( |4.1|) is just the statement that gij is group invariant: 

CQ^9ij = 0- (4.4) 
Using ( ^.1| ) together with the field equations 

] =0, (4.5) 



dL ( dL 



one may rewrite the invariance condition in the form 



(9/ dT 

= 6L = —6<p^ + = J.-,, 6C (4.6) 



where the ja'^ are the Noether currents: 



dL 

3-" = 7^Q'- = -^''9^J Q\ . (4.7) 



Expression ( |4.7| ) is valid for any sigma model. It is completely equivalent to the cor- 
responding expression obtained when the alternative action is used, involving the linear 
embedding variables (j)'^ and the Lagrange multipliers Xa- In this case the infinitesimal group 
transformation law takes the form 

6^ = G\t<p'SC (4.8) 

or, with Latin indices suppressed, 

5^ = G^^SC, (4.9) 

where the matrices Ga satisfy the commutation relation 

[G^,Gp] = G,c\p. (4.10) 

Since the Lagrange multipliers are invariant under the group they play no role in the defini- 
tion of the currents, and one finds 

dL 

Ja^ = Tr-Ga<P (4.11) 



where L is the Lagrangian for the alternative action. Expression ( [4.1 1| ) is usually more 
convenient to use than ( |4.7| ). 

Since the constants in eq. ([4 .61) are arbitrary we have 



Ja^,/. = 0, (4.12) 

which implies conservation of the total charges 

qa = j^Ja^dT.^. (4.13) 

Here S is any smooth deformation of a global spacelike Cauchy hypersurface in spacetime 
and dH^ is its surface element. 

In the case of the 0(l,2)/0(2) x Z2 model it is convenient to replace the single group 
index a by a pair of indices ah, so that equation ( ^4.8| ) takes the simple form 



5r = \G\db = 5i\ 0'' . (4.14) 

Here we have 

G\db = S\ Vdb - S\ Vcb, SCb = Vbc Sr ■ (4.15) 

The group parameters 6^,°'^ are those of an infinitesimal Lorentz transformation in the con- 
figuration space and satisfy 

S^'''' = -6^^\ (4.16) 
Reference to eqs. (|2.5| ) and ([4.11|) shows that the currents are 



Jab^ = -fi%d (p'^^Cabe = /i' {<Pa 06,^ - <Pb 0a,^) , (4.17) 

where (pa = Vab(p^- If the hypersurface S is chosen to be x° = constant then the corresponding 
charges can be expressed in the form 



qab 



J iab^ rf^X = -/i^ j (0„ 0b_o - 06 0a,o) . (4.18) 



5 Functional integrals. Ward-Slavnov identities. 

Since we shall ultimately attempt to define the quantum theory of the 0(l,2)/0(2) x Z2 
model by means of the functional integral of Feynman, it is appropriate that we discuss 
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quantum problems in functional- integral language from the outset. In principle the functional 
integral gives a direct representation of transition amplitudes: 



(out I in) = N j e^^I-^VM dcj) (5.1a) 
= N f e'^^''''^U(j)dX. (5.1b) 



Here |in) and |out) are vectors corresponding to arbitrary "in" and "out" states, the inte- 
grations are over field histories satisfying boundary conditions appropriate to those states, 
and it is understood that expressions ( p.l|) and ( |2.5| ) for S and S may need to be made more 
precise by imposition of integration limits as well as addition of boundary terms. and N 
are normalization constants and fi[(f)] is the configuration space measure: 



= Ua'^' i<t>ix)) , g = det (g,,) . (5.2) 

X 



In principle a measure functional A] should be inserted in the integrand of ( p.lb| ), but 
it would be simply a constant because the embedding space is fiat. The volume elements of 
( ^.la| ) and (|5.1b|) are respectively 

d(j) = lldcpXx), (5.3a) 

i,x 

d(f)dX = lllld(f)''{x)'[[dXA{x). (5.3b) 



Although expressions ( |5.2| ), ( p.3aD and ( [5.3b| ) are purely formal they have a ready interpre- 
tation on a lattice. 

The functional integral may also be used to represent "in-out" matrix elements of chrono- 
logically ordered operators: 



(out|T(A[0]) |in) = y A[0]e^^l'^V[0] # (5.4a) 
= N f A[<f)]e'^^'^'^U(j)dX. (5.4b) 



The symbol on the left of these equations is to be understood as a quantum operator, on 
the right as an integration variable. 

Suppose now the dummy variables of integration 0* in eq. ( ^.4a| ) are replaced by 

cj)'' = (l>' + Q\mC, (5.5) 
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where the 5^", instead of being constants as in eq. ( f4.1| ), are now scalar functions of compact 
support in spacetime. The replacement ( |5.5| ) leaves ( p.4a| ) unchanged, but it produces the 
following changes in the individual factors under the integral sign: 



Q\{4>{x))6C{x)d''x, 



(5.6) 



d X 



-%e 



(5.7) 



S {f,[cl>] dcjy) 



fJ'[(f)'] d(f)' — dcf) 



Y^g-'" [{dg^'Vdct^) Q\ + g^l''dQ\ld<\^\ 5C • (5.8 



The second line of ( ^.7|) follows from the observation that the Lagrangian would be invariant 
under ( ^.51) if the 5^" were constants. Expression ( |5.8|) actually vanishes by virtue of the 
following immediate corollary of eq. ( [4.41) : 

(5.9) 



= £q„ g'/' = d (g'/'Q\) /d<p^ . 



Taking note of the arbitrariness of the 5^" one may therefore infer 



= N 



\6(j)'{x 



-Q\ {<P{x))-zA[<P]f,^{x)\e^'^^ifi[<j,] 



(out| T ( ^^Q\ (0(x)) - ^A[^]f,,ix)] |in) . 



(5.10) 



\(50*(x) 

Because the "in" and "out" states are themselves arbitrary one can, in fact, infer the operator 
identity 

fl / A 4 f^l . \ 

(5.11) 



Note here that, because the functional integral of the difference of two functionals is the 
difference of their individual integrals, the chronological ordering that eq. (|5.4a| ) defines 
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commutes with the operation of differentiation with respect to the spacetime coordinates. 
Equation (|5.11| ) is a shghtly generahzed version of the standard Ward-Slavnov identity. 



The Ward-Slavnov identity can be extended to the embedding variables 0" by working 
with expression ( p.4b|) instead of ( |5.4a] ) and replacing ( |5.5| ) by 



0' = (1 + G„(5r)0. (5.12) 

The invariance of the volume element dcp under this transformation follows from 

trG« = 0, (5.13) 

which in turn follows from the semisimplicity of the invariance group. The extended Ward- 
Slavnov identity takes the form 

^ -T {A[<p]j^^{x)) = -iT { I^G, <P{x) \ . (5.14) 



dxf^ \6(j){x) 



6 Current algebra 



Setting = 1 in eq. ( |5.14|) one gets the operator conservation law 



?■ ^ =0 

On the other hand, choosing A[(f)\ = (f){x'), writing 

T {(j){x')j^'^ix)) = 9 - x°) </)(x')ja'^(x) + e (x° - ja^{x)<t){x') 

where 6 is the step function, and making use of ( p.lD , one gets 

'(p{x'),]J^{x)\5(x'^ -x^) = iG^(f){x)6{x' -x), 
which is equivalent to the equal-time commutator 

= 0(x)(5(x — x'). 



(^)) ja I 



Integration of this equation over x' yields 

[0, g„] = iG^ • 



(6.1) 



(6.2) 



(6.3) 



(6.4) 



(6.5) 
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which, by virtue of the group invariance of the Lagrangian L, in turn yields 



dL 



and hence 



■ dL 

-l-r—Gn 



^1 



Using (|4.15| ) one easily finds, for the 0(l,2)/0(2) x Z2 model, 

= 0, 



(6.6) 

(6.7) 
(6.8) 



'(f)\x),j,,\x')]6[x^ -x" 
'<P\x),jo.\x')]6{x'-x" 
<P\x),m\x')]6{x'-x" 



(6.9) 
(6.10) 
(6.11) 



(6.12) 



= i d'^i (f)j{x) - (j)i{x) 6{x - x') 

= i(j)i{x)6{x — x') , 

= —i6^i (j)o{x)6{x — x') 
= i6''i(f)%x)6{x-x'), 
where k G {1, 2}. 

It is sometimes convenient to use chart coordinates in configuration space other than 
the 0*. For the 0(1, 2)/0(2) x Z2 model these will be chosen so as to transform like the 0* 
under the 0(2) subgroup. The most general such coordinates are 

^' = f{\^\)<j>\ (6.13) 

where /(|0|) is a positive- valued smooth function such that {(pl f has everywhere 

positive slope. For example, if 



sinh ^ 101 
101 



then the are the Riemann normal coordinates a*: 

= s cos 9 
a"^ = s sin 6. 
It is readily verified that in general we have 

: i (fjix) — S^j ^i{x) 5{x — x') 



\p\x),3,f{x')\5(x'-x"') 
ip^{x),]Qi'^{x')\5(x^ -X' 



(6.14) 
(6.15) 

(6.16) 
(6.17) 
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7 Vacuum and 1-particle states 



The classical energy corresponding to the action ( p^.l] ) is nonnegative and vanishes for con- 
stant fields. Each different constant field corresponds to a different classical vacuum state. 
Since 0(1,2) acts transitively on the configuration space the classical vacua are obtainable 
from one another by group operations. One expects the degeneracy of the classical ground 
state to be refiected in a corresponding degeneracy of the quantum ground state. If |vac) is 
the vector corresponding to one of the quantum vacua then the vectors corresponding to the 
others are obtained through multiplication by unitary operators exp [iqa ^"). 

Since all points of configuration space are equivalent under 0(1, 2) one can always choose 
a Lorentz frame in the embedding space so that the classical vacuum state is 0^ = 0^ = 0. 
The corresponding quantum vacuum is fixed by the conditions 

01 = 02 ^ 0, (p" = (vac| 0" |vac) , a G {0, 1, 2} . (7.1) 

Here the normalization 

(vac I vac) = 1 (7.2) 

will be assumed. Each quantum vacuum will also be assumed to be stable. 

The degeneracy of the quantum ground state implies that the basic particles of the theory 
are two Goldstone bosons. The masslessness of these particles is already suggested by the 
form of the first term in the expansion ( |2.4| ) of the classical action. Since the classical 
action generates no bare vertices having fewer than four prongs, perturbation theory, if 
it were valid, would imply that these particles, despite being massless, are stable.0 Since 
perturbation theory fails one must either assume that these particles are stable or else assume 
that 2-point functions exist and that the quantum vacuum, like the classical vacuum, has an 
0(2) invariance that is never broken. In the latter case Goldstone's theorem will guarantee 
that these particles are stable, for, as will be shown later, the degeneracy of the quantum 
ground state is never removed. 

Since the S'-matrix of the theory is expected to be independent of the choice of inter- 
polating field one should in principle be able to construct particle creation and annihilation 
operators out of any of the fields (|6.13| ) including 0* itself. Asymptotic arguments imply 

^Phase space restrictions forbid the decay of a single massless particle into three or more others. The 
derivative coupling also suppresses the decay. 
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that the creation operators are given by 

a*(p)* = ifi^ J^^ u{x, p) - ^ j rfS^ , (7.3) 

where "+cxd" and "— oo" denote Cauchy hypersurfaces in the remote future and past respec- 
tively, either asymptotic region being usable because of the stability of the particles. Here 
u is a mode function that should strictly have the form of a wave packet but which, if the 
packet is big enough, may effectively be represented by a plane wave: 

u{x,p) = (27r)-3/2(^2cj)-i/2gip-x 

u = = \p\ , = 0. 

In the Appendix it is shown that if the 1-particle state vectors are defined by 

|i,p) = a* (p)* I vac) (7.5) 

then the normalization 

(^,P|J,P') = StjS{p-p') (7.6) 



(7.4) 



is secured by choosing the coefficient /x^ in (|7.3| ) to be the square root of the reciprocal of 
the residue at the "pole" at = in the Fourier transform of the 2-point function: 

(vac|T(,.(.)^^(.0) Ivac) = -^J -Jm^T^^ . (7.7) 

If were replaced by 0* and if all the terms in the integrand of ( p.4|) other than the first were 
missing, so that the theory were "free," then /i^ would be simply the bare Planck mass fi. 
For the nonlinear theory it must be computed. This, in fact, is a major goal of the numerical 
simulation. 

Strictly speaking, the singularity of the integrand of (|7.71 ) at = cannot be a pole 
when the massless particles interact but must be a branch point. The coefficient fi^"^ can 
nevertheless be identified because the unwritten terms ". . . " in the denominator of the 
integrand generally vanish more rapidly than as — 0. For example, in perturbation 
theory the first unwritten term behaves like p^lnp^. 

The proof in the Appendix, of the normalization ( |7.(j| ), makes use of the mode-function 
orthonormality relations 



-i u{x, p)* - ^ j u{x, p') dT.^ = 5(p - p') 

-I f u{x, P) - 1^ u{x, p') dE^ = 0, 



(7i 
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and of the defining condition 



a*(p) |vac) = 0, (vac| a*(p)* = 



(7.9) 



for the vacuum state vector. Equations ( [7.3|) and ([Tig ), together with the completeness 
of the mode functions, imply that (^"^{x) may, in the remote past or future, be effectively 
represented by 



(p\x) = / M(x,p)a*(p) + u(a;,p)*a*(p)* d^p, x^ ^ ±oo 



(7.10) 



This in turn implies 

(vac| if\x) |vac) = (7.11) 
for z° — » ±oo. Because the action is an even functional of the fields eq. (|7.11|) in fact holds 



for all X (cf. eq. (|7[l|)) 



8 Renormalization via the effective action 

Renormalization is most easily discussed in terms of the effective action, denoted here by 
or r A according as the classical action is taken to be S[(f)\ or S'[</), A]. Consider 
first f . The traditional way to define it is to introduce sources J a, J a coupled to the 0" and 
A^ respectively and to set 



iW[J] 



(out, vac I in, vac) = N U[s[4>A+S{Jar+JAX-)d^ ^^^^ 



where, because of the presence of the sources (which are assumed to have compact support in 
spacetime), one must now distinguish between "in" and "out" vacua. The functional W[J] 
is used to define the following quantities: 



Qai...a'^A'{...A'J' 



6 S 



6W 

JTa 

5 



-w ^ 



(8.2) 



.3) 



It is straightforward to show that these quantities are related to the so called n-point func 
tions: 



La jM 



iC' , etc. , 
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the n-point functions themselves being special cases of the general average 
{A[<p,X]) = 



(out, vac 


r(40,A]) 


in, vac) 


(out, vac 


in, vac) 



J ^i[Sl<t>,X]+ J {Ja r + JA A^) d-^x] 



(8.6) 



6" and A"^ are called mean fields and the G"s are known as correlation functions. 

The effective action is derived from the functional W[J] by the Legendre transformation 



0, a] = W[J] - I (.Ja r + Ja A^) d^x . 



It is not difficult to verify that F satisfies the equations 

6T 



-Ja. 



1.7) 



(8.8) 



Functional differentiation of these equations yields 



tG 



d x" = —6a S{x — x'), etc. 



(8.9) 



which reveals G"''', G"^', G"^^' , G"^^' as Green's functions of the second functional derivative 
of F. These are the full propagators of the theory. Further functional differentiation of 
eqs. ( |8.9|) enables one to express the higher order correlation functions in terms of these 
propagators together with the functional derivatives of F of order three and higher, which 
are known as the full vertex functions. The resulting relations have a tree graph structure 
that allows one to recognize F as the generator of the 1-particle irreducible amplitudes of 
the theory. 

The effective action F associated with the classical action S[(j)] is constructed in a 
different way. The reason for this is that although the variables 0" transform linearly under 
the invariance group the variables 0* transform nonlinearly among themselves. The direct 
coupling of the latter variables to sources would lead to a F which has no simple transforma- 
tion law under the group. What one does instead is first note that the above construction 
of F is completely equivalent to defining it implicitly (or recursively) by 



jru,A 



N / expil S[(j),X] + 



5T (f),X . . ST (f), X , ^ 



d X 



(j)dX. 
(8.10) 
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From this it is easy to see that F must be hnear in the A"' 



(l),X =T (j) + A^C, 



(8.11) 



One also notes that because the 0", A^ transform hnearly under the group so do the 0", A^: 

S^"" = G\b ^''SC, 6X^ = 0. (8.12) 

The sources J^, Jy^ may be assumed to transform contragradiently, and F and the Cj^ 
are therefore group invariant. 



The effective action F 



0, A in the same way as the classical 



is now constructed from F 
action S[(j)] is constructed from S[(f), A]. One sets the sources Ja equal to zero and solves the 
constraint equations 



5F 



0,A 



6X^ 







(8.13) 



for the superfluous embedding variables (or rather their barred forms) in terms of the chart 
variables in the coset space, and then substitutes into f to get F . In the case of the 
0(l,2)/0(2) X Z2 model, if the 0° are constant fields (e.g., when the sources Ja vanish) 
equation ( |8.13| ) is necessarily equivalent to 



(8.14) 



for some constant Z. 

If the model is to have a consistent continuum limit the effective action must, because of 
symmetry requirements, reduce at low energies to a local action of the classical form ( |2.1| ), 
but with possibly a new value of /i. Thus 



0] = -\-^' I {m) ^'J'^'d'x + AT [0' 



(8.15) 

where AF, which includes nonlocal contributions, becomes important at higher energies.^ 
The metric cjij is group invariant (i.e., it satisfies eq. ( [4.4|) ) but, because the 0" are constrained 



by (|8.14|) to a generally different hyperboloid than the original classical variables 0", Qij is 
not equal to the metric (|3.3|) but is instead given by 



gij — 5ij 



2 • 



(8.16) 



■^The clean separation of T into "leading" and "subleading" parts is implied by the same power-counting 
arguments that show the theory to be perturbatively nonrenormalizable. 
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The original metric is restored if the quantum operator fields are replaced by renormalized 
fields: 



" - Z-^cf)" . (8.17) 



The effective action then takes the form 

r = -i/i/j' / Qi, (0i?(x)) '^^^J'^^^ + Ar (8.18) 

where 

^^R = ZJ^. (8.19) 

/iij is the renormalized or "experimentally observed" Planck mass. Note that because the 
Lagrangian is nonpolynomial the replacement ( |8.17| ) renormalizes an infinite number of ver- 
tex functions at once. Note also that although ( p. IT) ) looks like a standard "wave function 
renormalization" the renormalization constant Z is not obtained from the residue of the 
pole of the 2-point function. Rather, it is the constant /i^ of eq. (17^71) that is defined by the 
residue. 

When the sources Ja vanish and the "in" and "out" vacua become identical, the vacuum- 
fixing conditions ( [7.1|) imply 

Z = <tP = (vac| / |vac) , (8.20) 

(vac| T (0ij(x)0^^(x )) |vac) = -iG% , (8.21) 

(cf. eq. ( ^.5| )). is the Green's function of the second functional derivative of the effective 
action with respect to the renormalized fields. The structure of the local part of expression 
( ^.18| ) implies that the low-energy behavior of this Green's function is given by 

G'^ = -±- ^ d^p. (8.22) 

^ (27r)4i ;,|(p2_,o) + ... ^ ^ > 



Comparison of eqs. (7.7), ( |8.21| ) and (|8.22|) allows one to infer 



/ii? = (8.23) 



Comparison of eq. (|8.19|) and (|8.23|) in turn informs one that the /i of eq. (^.15|) is just the 



of eq. ([77^ ) when = (f). 

9 Renormalization via current algebra 

The renormalization analysis above makes special use of the variables 0", which transform 
linearly under the group. Determination of the renormalized Planck mass /i/?, however. 
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should not depend on which field variables are used. The generality of eq. ( [7.7| ) suggests 
that the choice of variables should, in fact, be irrelevant and that for every /i^ there should 
be an easily computable such that 

^R = Z^Ji^. (9.1) 

In this section we show that this is indeed the case. 

Consider the vacuum average (vac| T (y9*(x')jaf,^(a;)) |vac) where is one of the fields 
( |6.13| ). Because of the displacement invariance of the theory this average must be a function 
of the differences x'^ — x'^ . It must also transform as a vector under Lorentz transformations. 
Its Fourier transform must therefore have the form 

"•(^'--■) (vac| T (ip\x')3ab^{x)) |vac) rf^x' = jfX\,{p') (9.2) 

where the X'^ab are certain scalar functions of p^. These functions may be determined by 
multiplying eq. ( |9.2|) by and invoking the Ward-Slavnov identity ( ^.14 ): 



p'X\, 



—I 



d 



dx'i^' 



^ip-{x'—x) 



' / e*^ (^'-^)-^(vac|T(y.^(x')ja6na;)) l^ac) rfV 



(vac 



T (ip\x')jab^{x)) |vac) (fx' 



dx^ 



I e*f("'-")(5(x' - x) (vacl ^G\m<P'^ |vac) rfV 



(vac I 



dip'' d^p^ 



ba I vac) . 



(9.3) 



The integration by parts and subsequent replacement of d/dx'^ by —d/dx^, leading to the 
second line, assumes that the vacuum average falls off sufficiently rapidly at infinity. Omission 
of the chronological ordering symbol in the last two lines is allowed by the ultralocality of 
the operators involved. 

X^ab is seen to have a simple 1/p^ behavior. The pole at = can be removed by taking 
the Laplacian of the original vacuum average and integrating by parts. We do this for the 
special case a = 0, b = j and note that, because the vacuum average depends only on the 
differences x'" — x'^, application of is equivalent to application of □: 

J ^^p■{x'-x)^, ^^g^^i J. |vac) 



-pV^'oj = (vac| ( - ) |vac) 



6ijp^ (vac I 



1 



/(I0l) + 2l'^l/'(l'^l) 



0° I vac) 



(9.4) 
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The last line is obtained by invoking the explicit form ( |6.13[ ) and the fact that the vacuum 
average must be proportional to 5ij since, as far as the indices i and j are concerned it is an 
0(2) invariant tensor. This follows from the fact that the vacuum is 0(2) invariant, being a 
zero-eigenvalue eigenstate of the charge qi2. 

Now consider the matrix element (vac| joj'^{x) \i, p). Since |vac) and \i, p) are normalized 
physical state vectors, and since jab'^ is the physical current that yields the physically observ- 
able charges Qab (which generate group transformations), this matrix element must depend 
only on observable renormalized quantities. Moreover, its value must be independent of what 
field variables one uses to calculate it. Making use of eqs. (|7.3| ), (|775|) , ( [7.9|) and (|9.4| ), one 
gets 



(vac| joj''(a;) =iiJ^ 



> \ I d d 



vac 



joj^{x)ip\x') |vac) dK 



oo J — oo 



u{x,p) I ^ ^ ^ I (^^^1 ^ ['P\x')kj^{x)) |vac) dT.'^ 



W - □') (vac| T U\x')]Qj^'{x)) |vac) rfV 



(vac I 



fm) + \\ci>\rm) 0ivac)%pM^,p). 



(9.5) 



In the special case / ( |0| ) = 1, y?' = 0*, //^ = //^ = /i, the final vacuum average is just Z 
(see eq. ( p.2(]| )), and it follows from eq. ( |8.23| ) that 



(vac| joj^(a;) |i,p) = -ifinSij p''u{x,p) 



(9.6) 



showing that this matrix element does indeed depend only on observable quantities. More- 
over, eq. (|9.1| ) now follows, with the identification 



= (vac| [/ ( 101 ) + I 101 /' ( 101 )J |vac) . (9.7) 
This identification allows eq. ( ^.21|) (with G% given by ( ^.22|) ) to be generalized: 

(vac|T(^y.)^,(.'))|vac> = J -J^^^ll^i',, (9.8) 

where 

'PR = Z-'^\ (9.9) 
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Note that since the fields 0^ and <y9^ are not identical the unwritten terms in the denominators 



of the integrands of expressions ( ^.221 ) and (|9.8| ) will not generally be the same. Only the 
coefficients of will be identical. 

We now have two completely independent ways of computing ^r- One is to compute 
by (|9.7|) , carry out the field renormalization ( |9.9|) , and determine the residue of the pole at 

= in the Fourier transform of the 2-point function ( |9.8|) . The other is to compute the 
matrix element ( |9.6|) . 

The first method can immediately be translated to Euclidean space and adapted for the 
computer. The second requires a little modification. Instead of working directly with the 
matrix element ( p.6| ) we arrange for it to appear in a sum over intermediate states: 



(vac|T(jo.nx)jo/(x'))|vac) 

= (vac| [9 (x° - J^^^{x)J,,\x') + e - x°) joi"(x')jo^^(x)] |vac) 

2 „ 

J2 - (vac| joi''{x) \k, p) {k, p| joj^ix') |vac) 



k=l ■ 



+ e - (vac| Joj'^ix') p) {k, p| joi^{x) |vac) ] d^p 



+ 



(9.10) 



Here we assume fi ^ u so that we may use the naive definition of chronological ordering 
without worrying about Schwinger terms. In the final sum the omission of the vacuum as an 
intermediate state follows from the Lorentz invariance and 0(2) invariance of the vacuum, 
which implies 

(vac| joi^(x) |vac) = 0. (9.11) 

The other unwritten terms in ( |9.10| ) are sums over intermediate states involving two or more 
particles. These are expected to be of higher order in than the written term when the 
Fourier transform is taken. 

Fourier transforms will here be defined by 



f{p) = (27r)-2 J f{x)e'P-^ d^x . (9.12) 
To obtain the Fourier transform of ( |9.10|) first insert (|9.6| ) into the sum, obtaining 
(vac I T {joi^{x)joj''{x')) I vac) 
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{2n 



2 e 



— 



which immediately yields 



(vac|T(Joi''(p)joi''(p')) |vac) = -in\5i. 



/i ^ I/, 



(9.13) 



/i^z/. (9.14) 



Since the unwritten terms are of higher order in it is evident that /i^ may in principle be 
determined by examining (Pl^ ) in the limit 0. 



10 Euclidean space 



To pass to Euclidean space one makes the variable replacements 



X 



-IX 



Po = 



(10.1) 



and carries out 90° rotations in the complex planes of these variables. In the functional 
integrals ( ^.la]) and ( |5.4a|) , as well as in expressions ( ^.2|) and (|5.3a|) for the functional 
volume element, the x'^ are simply labels, or generalized indices, on the dummy integration 
variables. Hence the replacements ( |10.1| ) would have no effect on the integrals themselves 
were it not for the fact that the integrands have an explicit dependence on x° through the 
dependence of the action functional ( |2.1|) on the volume element d'^x, which gets multiplied 



by —2. This dependence has the consequence that (|5.1a| ) and (|5.4a|) get replaced by 



(out I in) 



N 



-Se[ 



V[0] d(f), 



(out I A[(f)] I in) = N J A[0]e-^^1'^V[0] # , 
where Se is the Euclidean action, 

Se[(P] = Ifi'^ gij{(l){x))(l)\f,(f>'^^d^x, 



(10.2) 
(10.3) 



(10.4) 



and the Greek indices now run from 1 to 4. In eq. (|10.3| ) there may also be additional 
dependence on x^ and/or po if v4[0], |in) or |out) depends explicitly on one or both of these 
variables. In this case the corresponding dependence on x^ and/or p4 is obtained by analytic 
continuation. 

It will be noted in eq. ( p.0.4D that the Minkowski metric rj^i, has been replaced by the 
Euclidean 5^,^. Moreover, in eq. ( |10.3| ) the chronological ordering symbol T has been dropped 
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as irrelevant, the quantity on the left being now simply regarded as defined by the integral 
on the right. 

In Minkowski spacetime the Fourier components of |vac) behave like e^^^° in the 
remote past and those of (vac| (j)^[x) behave like e"**^^" in the remote future (see eq. (|7.10| )). 
The corresponding behaviors in Euclidean space are e^^^ and e~'^^* respectively, which imply 
that to get Euclidean vacuum averages one must integrate, in eqs. ( |10.2| ) and ( |10.3| ), over 
fields that vanish at infinity. On the computer, Euclidean space will be replaced by a 4-torus 
T^, the fields will satisfy periodic boundary conditions, and only their integrals over the 
4-torus will be required to vanish. The corresponding "Euclidean vacuum averages" will be 
indicated by the abbreviated notation. 

It is not possible to use a Euclideanized version of the action ( |2.5| ) on a computer because 
a straightforward Euclideanization of ( |2.5| ) leads to an action that is not bounded from below. 
One can use a Euclideanized version for theoretical purposes if the Lagrange multiplier field 
A is itself subjected to a rotation in the complex plane: 

\ = -i\E. (10.6) 

A Euclideanized effective action Ve (P,Xe can then be introduced, as well as a corresponding 



r 



E 



having the form 



gi, (0(x)) 0^',^ rf^x + ATe [4>] (10.7) 

(cf. eq. (8.15)). It will be left as an exercise for the reader to rewrite section 8 in Euclidean 
language. We remark only that the renormalization constants all remain unchanged, and 
equations (p.l7|) to ( |8.20|) take the forms 

rn = Z-'r, (10.8) 

= lfilJg.,{Mx))^'nJii^^d'x + ATE, (10.9) 

f^R = Zfx, (10.10) 

Z = 0° = (0°). (10.11) 



More generally, one may write 

fiR = Z^fi^ (10.12) 
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where 



/(I0l) + |l0lf (101 



and where /x^ may be identified from the relation 

1 



(2vr) 



r. . ip-(x-x') 
fiipP"^ -\ 



;io.i3) 



(10.14) 



which follows from ( [7.7|) and the fact that, in the transition from Minkowski spacetime to 
Euclidean space, the momentum volume element d'^p gets multiplied by z. In the special case 
in which the 99* are chosen to be the Riemann normal coordinates a* (see eqs. ( |6.14 ) and 
( |6.15| )) it is straightforward to verify that 



1 + 



s cosh s ' 
sinh s 



(10.15) 



The variables a* are particularly convenient to work with on the computer, and in the 
following paper an account will be given of the determination of fiR based on the equation 

6ij e*P-(^-^') 



1 



dp, 



a 



R 



a 



(10.16) 
(10.17) 



We note finally the Euclideanized versions of eqs. (|9.13|) and ( p.l4| ): 



f^R^ij f PjiPy_ ip.{x-x') li 



p^ 



d''p + 



{joMjojAp')) = f^R^ij'^^SiP + P') + ■ 



p^ 



^ z/. 



(10.18) 
(10.19) 



The Fourier transforms are here defined as in eq. ( |9.12|) , and the currents themselves are 
defined by 



Jab,,[X) = /i [(Pa(Pb,fj, - (Pb(Pa,t^j 



10.20) 



(cf. eq. i^)). 
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Appendix 

Equations ( [7. 3D and ( [7.5D imply 

(vac|(^*(a;) |j,p') = ifi^ J^^u{x' , p') (^-^) (vac| |vac) . (A.l) 
This may be rewritten in the form 

= -i/i^ / ) (□' - □') (vac| T {ip\x)ip^ {x')) |vac) rfV (A.2) 



since the insertion of the integral over the hypersurface +oo contributes nothing in view of 
( friOl) and the kinematics of the chronological product. Inserting expressions ( [7.4|) and ([7^ ) 
into (|A.2|) , and using the fact that u{x',p')D' = 0, one obtains 



(vac|v9'(x) |j,p') = (27r)-3/2^^(2cu')~'/' / c/V / d 



= (2.)-3/^;.,(2.')-/^/^^^%^^A 

= Sijfi^^u{yi,p') . (A.3) 

To reach the last line one makes use of the fact that the 5-function in the line above enforces 
the constraint p = p', = p'^ = 0, and hence the unwritten terms in the denominator of 
the integrand vanish. 

The adjoints of eqs. (|7.3|) and now yield 

{h P I J, P') = -ifiv J ^ u{x, P)* ( ^ - ^ j (vac| ip\x) \j, p') rfS^ 

= -.5., £^ p)* - ^ j u{x, pO d^, 

= S.,6{p-p'). (A.4) 
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